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蚊虫 驱 避 化 合 物 与 人 体 分 泌 的 引诱 物 
L- 乳 酸 分 子 间 弱 相 互 作用 的 计算 
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摘要 :【 目的 】 驱 避 剂 是 蚊虫 防治 的 重要 手段 之 一 。 前 期 研究 发 现 ,在 驱 避 剂 的 使 用 过 程 中 ,了 驱 避 化 
合 物 分 子 与 人 体 分 涛 的 引诱 物 分 子 存在 相互 作用 。 为 深入 了 解 这 种 相互 作用 的 本 质 ,本 研究 从 量 
子 化 学 的 角度 分 析 了 3 种 驱 避 化 合 物 分 子 与 引诱 物 分 子 工 -乳酸 相互 作用 过 程 中 几何 结构 、 结 合 
能 、 电 荷 布 居 和 热力 学 参数 等 量化 参数 的 变化 ,以 期 从 计算 机 模拟 角度 揭示 驱 避 化 合 物 分 子 与 引诱 
物 分 子 相互 作用 的 形式 、 强 度 和 规律 。【 方 法 】 本 研究 选择 了 3 种 在 驱 蚊 方面 具有 代表 性 的 驱 避 化 
合 物 , 即 8-ÉE X31 — $E ZEBE T SEES (RT) e P BE-1,2-99 — EAR BE( R2) de N,N-— G Ep v 
AA F BE (3 SUBE) ( DEET) ,使 用 Gaussian 09 量子 化 学 计算 程序 包 对 驱 避 化 合 物 分 子 单 体 RI1， 
R2, DEET 和 人 体 引 诱 物 分 子 单 体 工 -乳酸 以 及 驱 避 化 合 物 分 子 与 L- 乳 酸 的 缔 合 体 进行 结构 构建 ; 
在 M06-2X-D3/6-311G(d,p) 水平 上 优化 所 有 单 体 、 缔 合体 的 几何 结构 ,并 做 振动 频率 计算 ;在 M06- 
2X-D3/6-311 +G(d,p) 水 平 上 对 单 体 、 缔 合体 进行 单 点 能 计算 。[【 结果 】 驱 避 化 合 物 分 子 R1，R2 
fe DEET 与 上 -乳酸 的 缔 合 距离 在 2.64 ~2.74 A 之 间 , 缔 合 角 度 在 170 ~175° 之 间 , 缔 合 后 ,L- 乳 酸 
羧基 上 氧气 键 分 别 伸 长 0.019, 0.029 和 0.022 Å; RI, R2 fe DEET 与 -乳酸 缔 合 过 程 的 结合 能 
AE 4&7 - 50.278, — 72.385 和 一 66.977 kJ/mol, A X AH 2r 3| 7; —43. 715, -65.559 和 
—60. 465 kJ/mol,J8& X AS 分 别 为 —0.138, -0.159 fe - 0.157 kJ/mol - K, 自由 能 交 AG 分 别 为 
-2.652, -18.090 fe —13. 784 kJ/mol; R1, R2 和 DEET 5 L-ZLB& t Aat 42 P , L- $088 $9 ep 
由 0 e 分 别 变 为 -0.0248, -0.0250 4e -0.0168 e。[【 结 论 ) 计 算 结 果 表 明 , 了 驱 避 化 合 物 分 子 R1， 
R2 和 DEET 5 L- 和 乳酸 均 通 过 中 强 毛 键 缔 合 ,并 且 R2 与 L- 乳 酸 的 毛 键 作用 最 强 , 形 成 的 缔 合 体 也 
最 为 稳定 。 
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Calculation of the weak interaction between mosquito repellents and the 
human-secreted attractant L-lactic acid 
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Abstract: [ Aim] Using repellent is one of the most important methods for mosquito control. In our 
previous study, we found that there is an interaction between repellent molecules and human-secreted 
attractant molecules during the use of repellents. To understand the true nature of this interaction, in this 


study the interactions between three repellents and one attractant L-lactic acid were explored through 
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computational quantum chemistry based on the analysis of numerical changes of geometric structure, 
binding energy, charge population and thermodynamic parameters. This study might reveal the form, 
strength and law of the interactions between the repellent molecules and the human-secreted attractant 
molecules in silico. [ Methods 】 The structures of three representative mosquito repellents 8- 
hydroxylearveol formate ( R1), hydroxy citronellal-1, 2-propanediol acetal ( R2) and N, N-diethyl-3- 
methyl benzoyl amide ( DEET) , the attractant L-lactic acid, and their complexes were constructed by 
using Gaussian 09 software package. Geometry optimizations and vibration frequency calculations of the 
repellents, the attractant and their complexes were done at the M06-2X-D3/6-311G (d,p) level. The 
single point energy of the repellents, the attractant and their complexes were calculated at the M06-2X- 
D3/6-311 +G (d,p) level. [Results] The association distances and angles between the repellents R1, 
R2 and DEET and the attractant L-lactic acid were among 2. 64 -2.74 A and 170 - 175^ , respectively, 
and after the association, the bond length of oxygen hydrogen bond in carboxyl of L-lactic acid were 
stretched to 0. 019, 0. 029 and 0. 022 Å, respectively. The binding energies ( AE) of the association 
process were — 50. 278, — 72. 385 and — 66. 977 kJ/mol, the enthalpy change ( AH) - 43. 715, 
—65. 559 and -60.465 kJ/mol, the entropy change ( AS) 0.138, -0.159 and -0.157 kJ/mol - K, 
the free energy( AG) —2.652, — 18.090 and — 13. 784 kJ/mol, and the total net charge of L-lactic 
acid varied from 0 e to — 0.0248, -0.0250 and —0.0168e, respectively. [ Conclusion] The repellents 
R1, R2 and DEET and the attractant L-lactic acid are associated by moderately strong hydrogen bond. 
The hydrogen bond between R2 and L-lactic acid are the strongest, and the complex formed is the most 
stable. 


Key words: Mosquito repellent; L-lactic acid; weak interaction; hydrogen bond; binding energy; 


charge population; thermodynamic parameter 











蚊虫 可 传播 登革热 、 寨 卡 病 毒 病 和 西 尼罗河 热 
等 多 种 严重 疾病 (Bhatt et al., 2013; Nedjadi et al., 
2015; Petersen, 2015; Yakob and Walker, 2016) ,对 
AKEREKE K BS lU s 同时 随 着 全 球 气候 和 人 
类 生产 生活 方式 的 变化 ,蚊虫 防治 形势 日 趋 严峻 
( Hemingway, 2004; Murray et al., 2012; Lees et al., 
2014) 。 在 蚊虫 防治 过 程 中 , 驱 避 剂 发 挥 着 重要 作 
用 ( 陶 波 和 张大 伟 , 2014) ,新 型 驱 避 剂 的 研发 一 直 
是 研究 的 热点 。 

人 们 通过 呼吸 和 排 汗 释放 许多 对 蚊虫 具有 引诱 
活性 的 化 合 物 ,如 工 -乳酸 、 氮 1- 辛 炳 -3- 醇 等 (Acree 
ef al., 1968; Takken and Kline, 1989; Geier et al., 
1999) 。 蚊 虫 通 过 对 引诱 气味 的 溯源 ,可 以 精准 地 
找到 人 体 ,实施 吸血 行为 。 驱 避 剂 可 以 通过 不 同 的 
方式 阻止 蚊虫 对 人 体 定 位 和 叮咬 行为 ,包括 :(1) 阻 
断 蚊虫 对 引诱 气味 分 子 的 识别 ,导致 蚊虫 无 法 定位 
气味 来 源 ( Acree et al., 1968; Davis and Sokolove, 
1976 ; Ditzen et al., 2008) ; (2) 激活 蚊虫 对 驱 避 化 合 
物 的 感知 ,使 蚊虫 产生 厌恶 进而 主动 躲避 (Syed and 
Leal, 2008; Xia et al., 2008 ; Stanczyk et al., 2010) ; 
(3) 同 时 激活 蚊虫 多 个 嗅觉 重 白 ,使 蚊虫 产生 嗅觉 
紊乱 ( Bohbot and Dickens, 2010; Pellegrino et al., 
2011; Turner et al., 2011) 。 















































关于 驱 避 剂 对 蚊虫 嗅 沉 系统 影 响 的 研究 已 经 备 
受 关注 , 避 蚊 胺 CDEET) 涂抹 于 皮肤 上 与 不 和 皮肤 
接触 时 对 蚊虫 的 行为 影响 有 差异 ,而 且 DEET 与 L- 
乳酸 同时 存在 与 DEET 单独 存在 对 蚊虫 的 行为 影响 
也 有 差异 ( Dogan et al., 1999; Kline et al., 2003), 
说 明 DEET 与 人 体 分 泌 的 引诱 物 在 皮肤 表层 有 可 能 
存在 弱 相 互 作用 。 近 年 来 ,作者 所 在 课题 组 在 驱 避 
剂 及 其 作用 机 理 方 面 开 展 了 一 系列 的 研究 (Wang et 
al., 2008; 廖 圣 良 等 , 2012; Song et al., 2013; 忻 伟 
隆 等 , 2014; 许 锡 招 等 , 2014) , 提出 驱 避 化 合 物 与 
人 体 分 泌 的 引诱 物 可 以 发 生 相互 作 用 ,形成 缔 合 体 ， 
影响 蚊虫 的 叮咬 行为 。 但 是 其 具体 作用 形式 、 作 用 
强度 以 及 对 驱 避 剂 发 挥 驱 避 作 用 的 影响 有 待 进一步 
深入 研究 。 本 研究 拟 在 前 期 研究 的 基础 上 ,采用 量 
子 化 学 方法 从 几何 结构 参数 .结合 能 .电荷 布 居 以 及 
热力 学 参数 等 方面 分 析 驱 避 化 合 物 与 人 体 引 诱 物 分 
子 间 的 相互 作用 ,从 而 为 研究 新 型 蚊虫 驱 避 剂 及 探 
讨 蚊虫 驱 避 机 理 提供 理论 依据 。 


1 计算 方法 与 分 子 模型 


使 用 Gaussian 09 量子 化 学 计算 程序 包 ( Frisch 
et al., 2009) ,首先 对 驱 避 化 合 物 分 子 单 体 8- 羟 基 别 
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二 氧 欧 缕 醇 甲酸 酯 (R1) HEREN- ,2- 丙 二 醇 缩 
醛 (R2) DEET 和 人 体 引诱 物 分 子 单 体 [- 乳 酸 以 及 
驱 避 化 合 物 与 上 -乳酸 的 缔 合 体 进行 结构 构建 ; 然 
后 ,在 M06-2X-D3/6-311G (d,p) 水 平 上 对 所 有 单 
体 、 缔 合体 进行 几何 优化 和 振动 频率 计算 ,并 经 
Mulliken 电荷 布 居 分 析 对 缔 合 前 后 电荷 转移 状况 和 
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成 键 规律 进行 分 析 ; 最 后 ,在 M06-2X-D3/6311 +G 
(d,p) 水 平 上 对 所 有 单 体 、 缔 合体 进行 单 点 能 计算 ， 
进一步 计算 驱 避 化 合 物 分 子 与 引诱 物 分 子 缔 合 前 后 
的 各 热力 学 参数 的 变化 。R1，R2 ，DEET 和 工 -乳酸 
的 分 子 结构 以 及 R1，R2 和 DEET 与 L- 乳 酸 缔 合 体 
的 分 子 模型 如 图 1。 
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图 1 单 体 分 子 结构 与 缔 合 体 分 子 模 型 


Fig. 1 The structures of monomer molecule and complex 





A: 8- 产 基 别 二 氧 葛 缕 醇 甲 酸 酯 8-Hydroxylcarveol formate ( R1) ; B; 羟基 




















FIR- ,2- 丙 二 醇 缩 醛 Hydroxy citronellal-1 ,2-propanediol acetal ( R2) ; 


C; N,N- 二 乙 基 间 甲 基 茶 甲 酰胺 ( 驱 避 腕 )N,N-diethyl-3-methyl benzoyl amide (DEET) ; D: L-3L8£ L-lactic acid (LA) ; E; RI-LA 缔 合体 R1-LA 
complex; F; R2-LA 缔 合 体 R2-LA complex; G; DEET-LA 缔 合 体 DEET-LA complex. 








2 结果 


2.1 R1, R2 和 DEET 与 工 -乳酸 缔 合 
参数 

使 用 M062X-D3/6-311G (d,p) 方法 ,对 驱 避 
化 合 物 分 子 单 体 R1, R2, DEET 和 人 体 引 诱 物 分 子 
单 体 L- 乳 酸 , 以 及 驱 避 化 合 物 与 上 -乳酸 的 缔 合 体 进 
行 几 何 全 优化 和 振动 频率 计算 ,计算 结果 均 无 虚 频 ， 


的 几何 结构 


说 明 优 化 所 得 构 型 均 为 势能 面 上 能 量 最 低 点 ,是 最 
稳定 构 型 。 优 化 后 得 到 的 几何 结构 参数 列 于 表 1。 
L- 乳 酸 送 基 上 的 氧 氧 键 (0-H) 分 别 与 RI 和 DEET 
的 碳 氧 双 键 (C =0) 及 R2 缩 醛 环 上 的 碳 氧 单 键 (C- 
0) 缔 合 , 缔 合 距离 在 2.64 ~2.74 A 之 间 , 缔 合 角 度 
在 170 ~175° 之 间 , 根 据 对 毛 键 强 弱 的 一 般 性 判断 
( 吴 阳 等 , 2008; Arunan et al., 2011; 陈 静 , 2012), 
表明 R1，R2 和 DEET 与 工 -乳酸 通过 氢 键 缔 合 ,并 
且 缔 合 距离 在 2. 5 ~3. 2A 之 间 , 缔 合 角度 在 130 ~ 
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X1 RI, R2 和 DEET 与 工 -乳酸 缔 合 的 几何 结构 参数 


Table 1 Geometric parameters of the association between R1, R2 or DEET and L-lactic acid 











L-$LER323& Ef 0-H 键 长 (A) 
Bond length of O-H in carboxyl of L-lactic acid 











缔 合 位 点 缔 合 距离 (A) 缔 合 角度 (°) 
Association site Association distance Association angle 缔 合 前 缔 合 后 
Before association After association 
R1…LA -C=0:…H-0- 2.734 174. 335 0.966 0. 985 
R2--LA -C-O--H-0- 2. 641 170. 991 0. 966 0. 995 
DEET---LA -Cz0--H-0O- 2. 695 171. 284 0. 966 0. 988 








180*z7[H], JE -T FP vh ALES RUDI L-TLREXR AE E SUR 
键 键 长 为 0.966 À, 5 R1, R2, DEET 缔 合 后 , 键 长 
分 别 伸 长 为 0.985, 0.995 和 0.988 Å, 5 R2 缔 合 
时 , 伸 长 量 最 大 ,为 0.029 Å, RH R2 与 工 -乳酸 缔 合 
时 氧 键 最 强 ,DEET 次 之 ,RI1 相对 较 弱 。 
2.2 R1, R2,DEET 和 单 体 工 - 乳 酸 及 其 缔 合 体 的 
能 量 分 析 

驱 避 化 合 物 分 子 R1，R2, DEET 和 引诱 物 分 子 
L- 乳 酸 , 以 及 驱 避 化 合 物 与 上 -乳酸 缔 合 体 的 单 点 能 
() 列 于 表 2。 结 合 能 AE 即 为 缔 合 体 总 能 量 减 去 

















成 的 缔 合 体 最 稳定 ,RI 与 L- 乳 酸 缔 合体 的 稳定 性 
稍 弱 于 R2, DEET 与 L- 乳 酸 的 缔 合 体 , 这 和 几何 结 
构 分 析 结 果 一 致 。 
2.3 工 - 乳 酸 及 R1, R2 fü DEET 与 工 -乳酸 的 缔 合 
体 的 Mulliken 电荷 布 居 分 析 

表 3 列 出 了 了 驱 避 化 合 物 分 子 单 体 .引诱 物 分 子 
单 体 及 其 缔 合 体 的 电荷 分 布 状况 以 及 工 -乳酸 羧基 
上 氧 氢 键 的 成 键 电荷 布 居 。 由 于 工 -乳酸 单 体 呈 电 
中 性 ,所 以 它 的 总 净 电 荷 为 0, 当 工 -乳酸 与 驱 避 化 合 
物 分 子 R1，R2 和 DEET 分 别 缔 合 时 ,L- 乳 酸 的 总 净 





















































缔 合 前 驱 避 化 合 物 分 子 单 体 的 能 量 与 引诱 物 分 子 单 
体 的 能 量 之 和 。 计 算 公 式 为 :AE = Ee 7 ÈE repelen 
+ Ei )。 其 中 ,En 表示 驱 避 化 合 物 分 子 单 体 的 
量 ,Ei 表示 LL- 乳酸 的 能 量 ,表示 缔 合 体 的 
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R2 驱 避 化 合 物 分 子 单 体 R1, R2 和 DEET, 
引诱 分 子 单 体 工 -乳酸 及 其 缔 合 体 的 能 量 
Table 2 Energy of the repellents R1, R2 and DEET, 


the attractant L-lactic acid and their complexes 





E (kJ/mol) AE (kJ/mol) 
LA -902 069. 868 z 
RI -1 721 234.742 - 
R2 -1 933 927.494 z 
DEET -1 568 477.376 = 
R1:…LA -2 623 354. 888 -50.278 
R2--LA -2 836 069. 747 -72.385 
DEET---LA -2 470 614.220 -66.977 


R1, R2 fil DEET 与 L- 乳 酸 缮 合 的 结合 能 AE 
分 别 为 -50.278，-72.385 和 - 66. 977 kJ/mol, 一 
般 而 言 ,结合 能 AE 的 绝对 值 均 大 于 50 kJ/mol, 属 
于 强 氧 键 ( 陈 静 ,2012 ) ,表明 在 缔 合 过 程 上 能 量 是 
有 利 的 。 而 且 ,R2 与 L- 乳 酸 缔 合 的 结合 能 最 小 , 形 











电荷 由 0. 0000 e 分 别 变 为 -0.0248，- 0. 0250 和 
-0.0168 e, 说 明 在 缔 合 过 程 中 发 生 了 较 强 的 氧 键 
作用 ,而 且 R2 与 L- 乳 酸 分 子 间 和 毛 键 作用 强度 最 大 ， 
R1 和 DEET 相对 较 弱 ;电荷 转移 量 呈 现 这 样 的 规 
律 ,主要 是 由 于 RI 和 DEET 碳 氧 双 键 (C = 0) 的 存 
在 分 散 了 和 氧 原子 的 电荷 ,而 且 DEET 的 碳 氧 双 键 与 
葵 环 形成 共 斩 体 系 ,导致 电荷 转移 更 少 ,与 L- 乳 酸 
缔 合 的 氧 键 作 用 相对 较 弱 。 此 外 , 单 体 工 -乳酸 羧基 
上 和 氧 氧 键 成 键 电荷 布 居 为 -0. 0688 e, ^; R1, R2 和 
DEET 缔 合 后 ,成 键 电荷 布 居 分 别 为 - 0.0394, 
-0.0791 和 — 0. 0266 e, 表 明 工 -乳酸 与 R2 缔 合 后 ， 
L- 乳 酸 羧 基 上 氧 氢 键 减弱 ,与 RI, DEET HAJ, A 
氧 键 增强 ,所 以 工 乳 酸 与 R2 分 子 间 氧 键 作 用 最 强 ， 
与 电荷 转移 分 析 结 果 一 致 。 
2.4 缔 合 体形 成 过 程 中 的 热力 学 参数 

根据 振动 频率 计算 和 单 点 能 计算 得 到 驱 避 化 合 
物 分 子 单 体 与 L- 乳 酸 缔 合 过 程 的 部 分 热力 学 参数 
列 于 表 4。 驱 避 化 合 物 分 子 RI, R2 和 DEET 5j L- 
FURA ITEHIA ZE (AH) RIS AE CAS) 数据 表明 缔 
合 过 程 是 趋 于 有 序 的 放 热 过 程 。 自 由 能 变 (AG) 表 
HJ RI, R2 和 DEET 与 上 -乳酸 缔 合 过 程 可 以 自发 进 
行 ,自发 趋势 与 几何 结构 分 析 结 果 一 致 。 
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RI 引诱 物 分 子 单 体 工 -乳酸 及 驱 避 化 合 物 分 子 单 体 RI, R2, DEET 5 工 -乳酸 缔 合 体 的 净 电荷 布 居 和 成 键 电 荷 布 居 
Table3 Atom charges and Mulliken bonding population of the attractant L-lactic acid and the complexes of the 
repellents R1, R2 and DEET and L-lactic acid 








LA 羧基 O-H 成 键 电荷 布 居 (e) 


Bonding population of O-H in carboxyl of L-lactic acid 





LA 分 子 电荷 布 居 (e) 驱 避 化 合 物 分 子 电荷 布 居 (e) 
Charge population of L-lactic acid Charge population of repellents 
LA 0. 0000 - 
RI--LA — 0. 0248 0. 0248 
R2--LA — 0. 0250 0. 0250 
DEET---LA — 0. 0168 0. 0168 


表 4 缔 合 体形 成 过 程 的 热力 学 参数 
Table 4 Thermodynamic parameters 
in the formation of complexes 

AG (kJ/mol) 


AH (kJ/mol) | AS (kJ/mol - K) 





R1--LA - 2. 652 -43. 715 —0. 138 
R2--LA - 18. 090 - 65. 559 -0. 159 
DEET---LA - 13. 784 - 60. 465 - 0. 157 


3 讨论 与 结论 











驱 避 化 合 物 分 子 R1，R2 和 DEET 与 工 -乳酸 通 
过 氧 键 缔 合 ,并 且 均 属于 中 强 氧 键 。 缔 合 区 域 均 在 
驱 避 化 合 物 分 子 的 碳 氧 键 上 , 缔 合 距离 在 2.64 ~ 
2.74 A 之 间 , 缔 合 角 度 在 170 ~175° 之 间 , 并 且 缔 合 
后 工 -乳酸 羧基 上 氧 氢 键 依次 伸 长 0.019, 0.029 和 
0.022 Å; R1, R2 和 DEET 5 L- 乳 酸 缔 合 的 结合 能 
AE 依次 为 - 50.278, -72.385 和 - 66.977 kJ/ 
mol, AG 依次 为 -2.652，- 18.090 和 — 13.784 kJ/ 
mol ,表明 驱 避 化 合 物 分 子 与 工 -乳酸 缔 合 是 一 个 自 
发 的 过 程 ,并 且 形 成 的 缔 合 体 是 稳定 结构 ;Rl1，R2 
和 DEET 在 与 工 -乳酸 缔 合 时 ,电荷 由 驱 避 剂 分 子 向 
L- 乳 酸 转移 ,其 中 ,R2 与 上 -乳酸 缔 合 时 , 电 葵 转移 量 
最 大 ,由 于 三 者 的 分 子 结 构 和 原子 净 电 荷 分 布 不 同 ， 
在 与 上 -乳酸 缔 合 时 没有 表现 出 电荷 转移 量 递减 的 
规律 。DEET 与 L- 乳 酸 同 时 存在 与 DEET 单独 存在 
对 蚊虫 的 行为 影响 存在 差异 ,作者 所 在 课题 组 在 相 
关 研 究 中 发 现 , 莫 类 蚊虫 驱 避 化 合 物 RI 和 R2 也 有 
类 似 现 象 。 本 研究 计算 结果 表明 ,RI1，R2 和 DEET 
与 工 -乳酸 同时 存在 时 ,可 以 发 生 弱 相互 作用 , 且 表 
现 出 中 等 强度 的 氧 键 缔 合 ,与 驱 避 化 合 物 单独 存在 
相 比 ,在 几何 结构 参数 、 能 量 分 析 Mulliken 电荷 布 
居 和 热力 学 参数 方面 均 存 在 显著 差异 。 这 方面 差异 
与 对 蚊虫 行为 影响 差异 之 间 的 关系 值得 深入 研究 。 
本 研究 为 进一步 研究 人 体 皮 肤 上 了 驱 避 化 合 物 分 子 与 
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— 0. 0688 
— 0. 0394 
- 0. 0791 
— 0. 0266 


引诱 物 分 子 之 间 的 弱 相 互 作用 提供 了 前 期 研究 基 
础 ,同时 可 为 蚊虫 驱 避 作用 机 理 研 究 提 供 一 定 的 帮 
助 ,也 有 可 能 为 新 型 蚊虫 驱 避 剂 的 开发 提供 新 思路 。 
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